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ABSTRACT 

We have imaged the halo populations of a sample of nearby spiral galaxies using the Wide Field 
Planetary Camera 2 on broad the Hubble Space Telescope with the aim of studying the stellar population 
properties and relating them to those of the host galaxies. In four galaxies, the red-giant branch is 
sufficiently well populated to measure the magnitude of the tip of the red-giant branch (TRGB), a 
well-known distance indicator. Using both the Sobel edge-detection technique and maximum-likelihood 
analysis to measure the /-band magnitude of the red giant branch tip, we determine distances to four 
nearby galaxies: NGC 253, NGC 4244, and NGC 4945, NGC 4258. For the first three galaxies, the 
TRGB distance determined here is more direct, and likely to be more accurate, than previous distance 
estimates. In the case of NGC 4258, our TRGB distance is in good agreement with the the geometrical 
maser distance, supporting the the Large Magellanic Cloud distance modulus (m — M)q = 18.50 that is 
generally adopted in recent estimates of the Hubble constant. 

Subject headings: stars: luminosity function - stars: Population II - galaxies: distances - galaxies: 
individual (NGC 3031, NGC 253, NGC 4244, NGC 4945, NGC 4258) 



1. INTRODUCTION 

1.1. Spiral galaxy Halos 

The diffuse stellar halo component of galaxies represents 
a tiny fraction of the mass (about 1% in the case of the 
Milky Way; Morrison 1993) and the low surface brightness 
of the population makes study in extragalactic settings dif- 
ficult. However, galactic halos are unique laboratories for 
investigating fundamental galaxy properties. Age and/or 
metallicity distributions of halo stellar populations and 
their kinematics provide fossilized glimpses of the earli- 
est conditions of galaxy formation. Halo populations may 
give us answers to key questions about the chronology of 
galaxy formation: how the halo formation was related to 
the assembly of galactic mass and the formation of galax- 
ies. 

Different views have been advanced toward understand- 
ing the formation of galactic stellar halos. The first sce- 
nario was proposed by Eggen, Lynden-Bell, & Sandage 
(1962), where the authors proposed that the metal-poor 
stars in the Galactic halo were formed during a rapid col- 
lapse of a relatively uniform, isolated protogalactic cloud. 
This picture is generally viewed as advocating a rapid dis- 
sipative and monolithic collapse of protogalaxies. 

An increasing number of observational findings, both at 
high redshift and in nearby galaxies, together with cosmo- 
logical models of structure formation, challenge this gen- 
eral view, and suggest an alternative picture of how halos 
of galaxies may have formed. This alternative paradigm 
for halo formation is based on the idea that the galactic 
halo was formed via the accretion, over an extended pe- 
riod, of small metal-poor satellites which underwent inde- 
pendent chemical evolution before being accreted (Searle 



& Zinn 1978; see also Freeman 1987). The disruption of 
globular clusters, revealed by tidal tails (Grillmair et al. 
1995), may also contribute to the stellar halo (Aguilar et 
al. 1988). This is supported by the observation that halo 
field stars and globular clusters in the Milky Way have sim- 
ilar mean metallicities (Carney 1993). Note that this is not 
the case for M 31 stellar halo (Durrell et al. 2001, Perrett 
et al. 2002). An interesting consequence of this hypothe- 
sis is that a fraction of globular clusters may be stripped 
relics of highly nucleated dwarf satellites (see Freeman & 
Bland-Hawthorn 2002 for a detailed discussion). The dis- 
covery of the Sgr dwarf galaxy (Ibata et al. 1994), streams 
and moving groups in halos of both the Galaxy (Yanny et 
al. 2000; Ivezic et al. 2000; Dohm-Palmer et al. 2001) 
and M31 (Ibata et al. 2001) encourage renewed interest in 
halo formation through accretion of smaller protogalactic 
structures. 

Although great effort has been dedicated in recent years 
to investigate stellar populations and/or kinematic of Lo- 
cal Group spiral halos, it is unlikely that the problem 
of galactic halo formation may be understood entirely by 
studies of these galaxies alone. Questions about the uni- 
versality of halo formation mechanism(s), the similarity of 
halo global properties and their sensitivity to disk and/or 
bulge properties, among others, will require a larger sam- 
ple of galaxies. 

One means of addressing these questions is to extend 
the study of galaxy halos to other spiral galaxies. Such 
studies are now within reach of the HST. Our purpose is 
to study systematically the global properties of population 
II halos of nearby normal spiral galaxies. Such a program 
may reveal whether there was a universal formation mech- 
anism, or whether there are mixture of different stellar 
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halos; if the halos are dominated by formation of the nu- 
clear bulge and perhaps enriched by gas outflow, or com- 
posed of debris of disintegrated dwarf galaxies revealed by 
the presence of a large fraction of intermediate age stellar 
populations. Also of great interest is the investigation of 
whether halo and disk/bulge properties are correlated. 

In order to interpret the color-magnitude diagrams 
(CMDs) of resolved stars in nearby galaxies it is important 
to have reliable distance estimates. Fortunately, some- 
times the CMDs themselves provide such an estimate. 

1.2. A distance indicator to nearby galaxies: the tip of 
the red giant branch 

Stars on the first-ascent red giant branch (RGB) climb 
this phase with an expanding convective envelope and an 
hydrogen burning shell. While evolving through the RGB, 
low-mass stars develop an electron-degenerate core, which 
causes an explosive start of the core-He burning phase, 
the so-called He-flash, almost independently from the ini- 
tial stellar parameters, such as initial mass or abundance 
(Chiosi et al. 1992). The He-flash is followed by a sudden 
decrease of the stellar luminosity due to the extinction of 
the H-burning shell. Hence low-mass stars will accumu- 
late along the RGB, reaching their maximum luminosity, 
during this phase, at the tip of RGB (TRGB). This behav- 
ior translates to an abrupt discontinuity in the luminosity 
function (Rcnzini 1992). As the evolution of the RGB 
stars depends essentially on the core-mass, the number of 
stars per luminosity interval will be related directly to the 
He-core growth rate (Iben & Renzini 1983). Theoretical 
and observational investigations find the luminosity func- 
tion of RGB stars to follow a simple power-law (sec Zoccali 
& Piotto 2000 and references therein). The TRGB mag- 
nitude is the location where the RGB luminosity function 
truncates. Below metallicities [Fe/H] f« —0.7, the location 
of the TRGB in the / band is expected to be quite in- 
sensitive to age and/or metallicity (e.g.. Salaris & Cassisi 
1997). At higher metallicities, bolometric corrections in 
the I band become important and the TRGB magnitude 
becomes increasingly sensitive to metallicity (e.g., Salaris 
& Girardi 2005). 

The TRGB method is also well supported observation- 
ally. The location of the absolute I-band TRGB lumi- 
nosity for a sample of Galactic globular clusters spanning 
a large range of abundances, —2.1 < [Fe/H] < —0.7, over 
ages spanning 2-15 Gyr (Da Costa & Armandroff 1990), 
is quite stable and insensitive to age and metallicity (for 
[Fe/H] < —0.7), changing by less than ~ 0.1 mag (Lee et 
al. 1993). 

Good agreement is generally found between the dis- 
tances obtained using the TRGB method and classical dis- 
tance indicators such as Cepheid and RR Lyrae variables 
(e.g., Sakai et al. 1996, Ferrarese et al. 2000). Being a pop- 
ulation H distance indicator, unlike the Cepheid method 
(population I distance indicator), the TRGB method can 
be applied to any class of galaxies, making this an attrac- 
tive and easy technique to use to estimate galaxy distances. 

Because the outer regions of galaxies are devoid of signif- 
icant concentrations of gas and dust and the stellar density 
is low, studies of standard candles in the halo provide a 
means of estimating distances free of assumptions concern- 
ing internal reddening in galaxies or crowding. In § 2, the 



details of observations using the HST Wide Field Plane- 
tary Camera 2 (WFPC2) and data reduction are reported. 
In § 3, we present the detection of RGB stars and the meth- 
ods for measuring the TRGB using the 7-band luminosity 
function that we have employed. Finally, in § 4, we present 
the color-magnitude diagrams, luminosity functions, and 
derived distances for the galaxies studied. 

2. DATA 

Until recently, no spiral galaxy halo outside of the Lo- 
cal Group had been imaged to a depth sufficient to per- 
mit study of the halo metallicity distribution function or 
other stellar population studies, such as the presence of 
intermediate-age stars. HST observations have begun to 
change this. We have observed a sample of spiral galax- 
ies with the aim of unambiguously resolving their stellar 
halos down to one or two magnitudes below the TRGB. 
Our selection criteria of the sample galaxies are as follows. 
We seek spiral galaxies with (i) morphological types be- 
tween Sa to Sc, (ii) distance moduli (m — M)o < 29.5, 
(iii) high inclinations, i.e., i > 45 degrees, to reduce any 
contamination from the outer disk stars, and (iv) absolute 
magnitudes Mb < —18. The resulting sample of eight 
galaxies represents most of the nearby luminous inclined 
spiral galaxies. In this paper, we estimate the distances 
to four galaxies in our sample, i.e., NGC 253, NGC 4244, 
NGC 4258, and NGC 4945, using the tip of the red giant 
branch technique. Fig. 1 shows the locations of the ob- 
served halo fields, superimposed on the Digitized Sky Sur- 
vey images of these galaxies. Madore & Freedman (1995) 
argue that in order to minimize potential biases and to se- 
curely identify the TRGB, there should be at least 50- 100 
stars in the first magnitude interval of the RGB star lumi- 
nosity function. With sparser data sets the distance modu- 
lus tends to be overestimated. Four galaxies in our sample 
have luminosity functions that are too poorly populated 
to yield rehable TRGB distances: NGC 55, NGC 247, 
NGC 300, and NGC 3031. Hence these galaxies are not 
discussed in the present paper; the color-magnitude di- 
agrams of the observed halo fields for these galaxies are 
shown and discussed in paper III of this series. Tikhonov, 
Galazutdinova & Drozdovsky (2005) estimate TRGB dis- 
tances for three of these galaxies using larger HST data 
sets. 

The galaxy sample is summarized in Table 1. The galax- 
ies were observed through the F814W and F606W filters. 
Exposure times through the F814W filter were set to reach 
S/N = 5 in the WF camera for an absolute magnitude 
Mi = -1 (for galaxies with (m - M)o < 27) or Mi = -2 
(for NGC4258 and NGC4945, which are more distant). 
The F606W exposure times were set to reach the same 
S/N at the same absolute magnitude for metal-poor RGB 
stars. The exposures were typically dithered over 3 sepa- 
rate pointings (the range being from 2 to 6) with a pattern 
extending over approximately 0.25 arcscconds to allow re- 
jection of hot pixels and detector artifacts. 

The images were reduced through the standard HST 
pipeline, using the latest fiatfield observations and us- 
ing contemporaneous super-dark reference frames. The 
dithered frames were combined with iterative cosmic ray 
rejection using software in the stsdas. dither package based 
on the drizzle algorithm of Fruchter and Hook (2002). 
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Briefly, the images were drizzled (shifted, geometrically 
corrected, and resampled) to a common frame to construct 
a median image reasonably free of cosmic rays. This image 
was geometrically transformed back to the pixel grids of 
the original image and used as the truth image for identi- 
fying cosmic rays. The images were then re-drizzled onto a 
pixel scale of 0.1", masking out the cosmic rays, hot pixels, 
and detector artifacts, using optimal weights based on the 
counts at the sky level in each image. A value pixfrac=0.8 
was used to minimize the degradation of the image quality 
due to the resampling. 

The stellar magnitudes were measured through circular 
apertures with a radius of 0.15". These aperture mag- 
nitudes were corrected to total magnitudes using Tiny- 
Tim model (Krist 2004) point-spread functions. Aperture- 
corrections and zeropoints are applied separately for each 
CCD chip. Magnitudes are in the Vega system (i.e., Vega 
has magnitudc=0 through both filters, when measured 
through a 3"aperture) . Diffuse sources have been excluded 
from the catalogs based on the magnitude difTercncc be- 
tween measurements through 0.09"and 0.3"apertures. Ob- 
jects where this difference is greater than 1.3 magnitudes 
arc generally galaxies or blends of stars. This demarca- 
tion was determined by starting with a Tiny-Tim predic- 
tion of the ratio of fluxes for a point source measured in 
those apertures. We mark objects in the images whose 
flux growth exceeds that prediction and iterate to arrive 
at a value which preserves point sources while excluding 
those with extended radial profiles. 

The stellar photometry was corrected for charge-transfer 
efficiency (CTE hereafter) effects using the 31 May 2002 
version of the Dolphin (2000) equations. The CTE correc- 
tion amounts to a maximum of ^ 0.2 mag for faint stars at 
the top of the chip in fields with the highest background, 
such as NGC3031, while reaching ~ 0.5 mag in fields with 
low sky background, such as NGC55. Typically the max- 
imum correction for the data is —0.5 mag. 

Transformation of the instrumental magnitudes to stan- 
dard V and I magnitudes followed the prescriptions of 
Holtzman et al. (1995). First, we have corrected from the 
foreground extinction using Agoew =■ 2.677 x E(B — V) 
and A814W = 1.815 x E(B — V), then given the instrumen- 
tal magnitudes, we have determined the standard Johnson- 
Cousins magnitudes by solving iteratively the second de- 
gree equations relating the WFPC2-to-VI magnitudes us- 
ing the coefficients tabulated by Holtzman ct al. (1995). 
Fig. 2 shows the evolution of the resulting photometric er- 
rors for I-band as a function of the apparent magnitude. 
The typical trend of increasing photometric errors with 
apparent magnitude can be seen, where most of the stars 
observed in the fields are affected by photometric errors of 
a few tenths of a magnitude. 

3. MEASURING THE TRGB 

To avoid stellar populations that may bias the analysis 
of the luminosity function, we have restricted our analy- 
sis to stars that match the location of the RGB sequence 
of metal-poor stellar populations. To do so, we have re- 
tained only stars that are bluer than a certain cut. For all 
the galaxies discussed in this paper, the remaining num- 
ber of red giant stars after the color cut is sufficient for 
a reliable TRGB estimate by the criterion of Madore & 



Freedman (1995); i.e., there arc more than 100 stars in 
the one-magnitude interval fainter than the TRGB. It is 
worth mentioning that as long as the aim is measurement 
of the TRGB, there is no need to include all the observed 
stars to construct the luminosity function. When the red 
giant branch is not well populated, metal-rich stars, i.e., 
[Fe/H] > —0.5, that occupy the red part of the CMD where 
no calibration is available, might make the tip of the red 
giant branch less well identifiable. Our simple color cut 
accomplishes this goal, but will lead to different overall 
luminosity functions for different metallicity distributions. 
This is not a problem, provided that the TRGB feature 
itself is clearly visible. The mean metallicities of halo stel- 
lar populations of all our sample galaxies are at or below 
[Fe/H] < —0.6 (see the second paper of this series for more 
details), so any contamination from metal-rich stars is ex- 
pected to be small. 

We used two different methods to detect and to measure 
the TRGB I-band magnitude from the luminosity function 
of a stellar population. In the following subsections we 
present briefly our methodology. 

3.1. Edge Detection 

The most widely used procedure to estimate the TRGB 
I-band magnitude uses the abrupt termination of the RGB 
luminosity function at the TRGB magnitude. The TRGB 
discontinuity causes a peak in the first derivative of the 
observed stellar magnitude distribution (Lee et al. 1993; 
Madore & Freedman 1995). Sakai et al. (1996) update 
this method by Gaussian-smoothing the luminosity func- 
tion to avoid the arbitrary choice of the binning size, and 
use a continuous edge detection function. 

Due to the power-law form of the RGB luminosity func- 
tion, we use a kernel-smoothed logarithmic luminosity dis- 
tribution function filtered by a Sobel kernel to determine 
the luminosity of the TRGB; the edge detection function 
is: 

E{m) = logio(«'(TO + am)) - logio($(TO - a^)). (1) 

where am is the photometric error, estimated as the 
mean of photometric errors shown in Fig. 2, within a bin 
of ±0.05 magnitude about a magnitude m (Sakai et al. 
1996). The edge detection filter essentially measures the 
slope of the luminosity function, and hence it is very sen- 
sitive to the noise; any sudden change in the signal is ex- 
tremely amplified by the filtering scheme. To suppress any 
detection of statistically non-significant peaks in the first 
derivative of the luminosity function, we apply a weighting 
scheme to the edge detection filtering output (Sakai et al. 
1997; Mendez et al. 2002). Each E{m) is weighted by the 
Poisson noise (i.e. ,$(771)^/^). 

Cioni et al. (2000) have found the edge detection tech- 
nique to be a biased estimator of the TRGB magnitude. 
However, the bias depends on the photometric errors and 
the amount of smoothing in the luminosity function at the 
location of the TRGB. For our data set, the bias is less 
than 0.03 - 0.04 magnitude (Cioni et al. 2000), much 
smaller than the random errors; therefore the distance 
modulus to the galaxy sample we derive using the edge 
detection technique would not be significantly affected in 
a comparison with the maximum likelihood (see below). 
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To estimate the statistical uncertainty in measuring the 
/-band TRGB magnitude, we use bootstrap resamphng 
of the data to generate a large number of samples drawn 
from the original data set to simulate the act of observing 
multiple times. Using this procedure to quantify the for- 
mal errors in measuring the TRGB has the advantage of 
using the empirical distribution function derived from the 
data itself, rather than a particular functional form (e.g., 
Gaussian) for the photometric errors. 

For each galaxy in our sample, wo perform the boot- 
strap resampling and calculate the statistic of interest 
N 1^ n X [log(n)]^ times, where n is the number of ob- 
served stars within the magnitude range where the lumi- 
nosity function model is valid in the case of maximum 
likelihood analysis, and the number of stars used to con- 
struct the luminosity function in the case of Sobel edge 
detection analysis. To this extent the bootstrap version 
of a statistic's sampling distribution matches the asymp- 
totic sampling distribution (sec Babu & Feigelson 1996 
and Efron & Tibshirani 1986 for reviews of the theory and 
application of this technique). In these bootstrap samples, 
the standard deviation of the I-band TRGB distribution 
is taken as the random error in the measurement of the 
TRGB magnitude. This procedure allows a more robust 
estimate of the statistical errors than using the full- width 
at half maximum of the observed peak in the first luminos- 
ity function derivative, as usually done (Sakai et al. 1996; 
see also Mendez et al. 2002 and Cioni et al. 2001). 

3.2. Maximum Likelihood Analysis 

An alternative to filtering is to fit models to the data 
using the technique of maximum likelihood; this avoids 
concerns over binning and kernel smoothing. While theory 
does not provide strong guidance on the exact functional 
form of the luminosity function near the TRGB, a simple 
two-power-law model appears to provide a good empirical 
fit over a restricted range of magnitudes near the TRGB 
itself. The TRGB itself is identified as the break between 
the two power laws. 

In this technique, the underlying luminosity function is 
treated as a probability distribution function (PDF) and 
we compare relative likelihoods of drawing the observed 
data set from the model PDFs as we vary the model pa- 
rameters. An advantage of this procedure is that no data 
binning is needed. In practice we have binned the both 
the probability distributions and the data in intervals of 
0.002 mag, as it simplifies the fitting procedure. 

Given a model of the stellar distribution on the RGB, 
$(m|^)k, where 6 is the set of free parameters characteriz- 
ing the model, the likelihood that the observed data set is 
drawn from this model is given simply by the product of 
individual probabilities that star i has an observed mag- 
nitude rrii within a range dm,: 



(2) 



where Pi is the individual probability defined as the ra- 
tio of the number of stars actually observed at a given 
magnitude by the total number of stars expected given a 
luminosity function model ^{rn\9). 

Photometric errors will affect the magnitude distribu- 
tion of the observed stars. To account for this, we convolve 



the intrinsic luminosity function models with an appropri- 
ate broadening function, S(to), to obtain the distribution 
of observed magnitudes. This function describes the prob- 
ability that a star with an intrinsic magnitude mint is ob- 
served to have the magnitude mots = "niinL + Thus the 
distribution of the observed magnitudes is: 

$(mj) = j g{m\6)E{mi,m,)dm (3) 

where g{m\9) is the intrinsic luminosity function model. 
We represent the measurement error by a Gaussian func- 
tion with a total dispersion cjm- That is, we have 
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where is the photometric error at the magnitude 
mi. The probability of observing a star of a magnitude 
mi in a magnitude range dm is simply g{mi\6) dm. The 
photometric errors smear the probability to the shape of 
the error function, centered at mt; then the probability be- 
comes E(mi, m) g{mi\9) dm. The photometric errors were 
determined from artificial star tests. We take the output 
of these simulations in 0.5 mag intervals for input stars of 
a fixed color V-I = 1.0. In each 0.5 mag interval, we de- 
termine the mean and standard deviation of the residuals 
of recovered minus input I-band magnitudes. For smear- 
ing the model luminosity functions, we use a set of Gaus- 
sian kernels with the mean and standard deviation that 
were determined from the simulations and with an overall 
normalization set by the fractional completeness in that 

0. 5 mag interval. The model luminosity fimctions are con- 
volved with these kernels prior to estimating the logarithm 
of the likelihood. 

The maximum likelihood fitting is restricted to the mag- 
nitude range where the effects of both the photometric 
errors and the incompleteness are not large. To avoid con- 
tamination from metal-rich stars that may bias the esti- 
mate of the TRGB location, we consider only blue stars, 

1. e., stars within the color range —2 < (V — I) < 2 (in prac- 
tice there are few stars bluer than V — I = 0.5). Ideally, 
it would be best to fit a curved TRGB to a locus of RGB 
models of all metallicities. However, in practice the reso- 
lution of existing isochrones near the TRGB is insufficient 
for this purpose, and also the isochrones themselves have 
been constructed via a series of interpolations that become 
increasingly uncertain near the TRGB. For this reason, 
the TRGB distance indicator remains an empirical test: 
a search for a clear jump or inflection in the luminosity 
function. 

To approximate the intrinsic luminosity function 
g{m\9), we use a broken power-law model (Cioni et al. 
2000; Mendez et al. 2002; Gregg et al. 2004). The free 
parameters characterizing the intrinsic model are the loga- 
rithm of the amplitude of the discontinuity at the TRGB, 
the apparent magnitude of the TRGB break, and the 
slopes of luminosity function on either side of the break. 

We used both the Ncldcr-Mcad (1965) simplex algo- 
rithm and "simulated annealing" (e.g.. Metropolis, 1953; 
Press et al. 1992) to determine the best-fit model. Re- 
sults of the two methods were similar, and well within the 
uncertainties derived from the subsequent error analysis. 
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To determine the uncertainties in the TRGB magnitude, 
we have stochastically sampled the full model parameter 
space by building a Monte-Carlo Markov Chain (MCMC) 
using the Metropolis-Hastings update mechanism. There 
is a growing body of literature on this technique (e.g., 
Knox, Christensen & Skordkis 2001; Verde et al. 2003), 
which provides an efficient way to characterize the likeli- 
hood manifolds in a multi-dimensional parameter space. 
The Metropolis-Hastings MCMC algorithm is a restricted 
random walk through parameter space. At each step in the 
random walk, magnitude and direction of the next step is 
drawn from a normalized "proposal density" q{y, x), where 
the new state y depends only on the current state x. The 
candidate y is accepted or rejected with a probability 

. f-. p{y) q{y,x)\ 

mm 1,-— — , (5) 

V pyx) lyx, y) J 

where p{y)/p{x) is the likelihood ratio for state x relative 
to state y. After a period of "burn-in," the distribution 
of locations in parameter space sampled by the chain con- 
verges to the "posterior distribution" p. Thus likelihood 
contours can be constructed simply by drawing contours 
of the density of points sampled by the chain. 

For our proposal distribution, we have used a Gaussian 
with (T = 0.1 in all of the fit parameters. We have also re- 
stricted the chains to operate within a "reasonable" range 
of parameter space: break amplitudes within a factor of 
100 of the best fit, TRGB magnitudes within 1.5-2 mag 
of the best fit, and LF slopes between —5 and -1-0.2. The 
chains were burned in for 1000 steps, and then 10000 steps 
were used to sample the posterior distributions. The prob- 
ability distributions thus derived are shown in the lower- 
right panels of Figures 3-6. We have verified that the re- 
sults are insensitive to small changes in the range of mag- 
nitudes being fit and the parameters used to define the 
MCMC. The confidence intervals (and indeed the best fit 
values) are nevertheless subject to a number of implicit 
prior assumptions. These include the assumption that the 
distribution can be adequately described by two power 
laws with a break, over the full magnitude range being 
fit. In practice, we do see small but significant differences 
in the confidence intervals as we change the magnitude 
limits and color limits. However shifts of 0.2 mag in these 
limits in all of our tests shift the peak of the likelihood 
distribution by less than the la widths show in Figs. 3-7. 

4. TRGB DISTANCES 

As a first step to investigate the properties of field 
stellar populations in spiral galaxy halos, we have deter- 
mined distances to a sample of galaxies. A summary of 
TRGB distances and uncertainties is presented in Table 
1. Color-magnitude diagrams, luminosity functions, and 
other properties for our galaxy sample are presented in 
the next subsections. The best fitting luminosity func- 
tion model, $(m|^), where 9 is the free parameter set that 
maximize the likelihood, is overplotted on the observed lu- 
minosity function for each galaxy. Note that the smooth- 
ing by the photometric errors, accounted for in the maxi- 
mum likelihood analysis, converts the abrupt break at the 
TRGB magnitude into a continuous RGB star distribution 
near the TRGB. 

Once the I-band magnitude of the TRGB is measured, 
we use the semi-empirical calibration given by Lee et al. 



(1993) to determine the distance modulus, starting from 
the relation: 

(m - M)/,o = Itrgb - Mhoi,TRGB + BCi (6) 

The bolometric magnitude, Mhoi,TRGB, and the bolo- 
metric correction, BC/, are function of the metallicity and 
the color of RGB stars (Da Costa & Armandroff 1990): 

Mboi,TRGB = -0.19 X [Fe/H] - 3.81 (7) 
SC/ = 0.881 -0.243 X (V-/)TKGi3 (8) 

where {V — I)trgb is the color of stars at the TRGB 
magnitude. The metallicity [Fe/H] is function of the RGB 
star colors (Lee et al. 1993): 

[Fe/H] = -12.65 + 12.6 x (^ - /)_3 5 - 3.3 x (F - /)!.3 5 

(9) 

where (V — /)-3.5 is measured at an absolute I-band 
magnitude of —3.5 mag. To measure {V — I)trgb^ we cal- 
culate the color distribution of stars within Itrgb ±0.1, 
and fit a Gaussian to the distribution. A similar procedure 
was repeated for the simulated bootstrap stars to measure 
the uncertainties, and a similar procedure was adopted to 
measure iy — /)-3.5. The absolute I-band magnitudes are 
given for each galaxy in Table 1. 

We assume that internal extinction is negligible as we 
are dealing with halo regions that are most likely dust free. 
To correct for the foreground extinction, we use the lOO/itm 
DIRBE/IRAS all-sky map of Schlcgel et al. (1998). 

4.1. NGC 253 

NGC253 is a well studied nearly edge-on Sc-type L* 
galaxy in the Sculptor group; it is an archetypal starburst 
galaxy, and is the closest edge-on infrared luminous galaxy 
(Radovich et al. 2001). Extended soft X-ray emission 
from a hot gaseous component was detected in the halo 
of NGC 253. This X-ray emission seems to be powered 
by a large-scale outflow, driven by the nuclear starburst 
activity (Pietsch et al. 2000; Strickland et al. 2002). 

The left panel of Fig. 3 displays the color-magnitude di- 
agram of the observed halo field. The stars in the halo 
are predominantly red; the main feature of the color- 
magnitude diagram is the RGB structure, which finishes in 
a well-defined TRGB. The red giant branch is prominent 
and wide, indicating a large spread in the metallicity dis- 
tribution of the halo stellar population. The stars above 
the TRGB and extending up to F — / ^ 4 may be Asymp- 
totic Giant Branch (AGB) stars covering a range of ages 
and metallicities. 

The color-magnitude diagram shows no indication of the 
presence of blue early-type main sequence stars in the ob- 
served halo field of NGC 253, as was reported previously 
by Comeron et al. (2001). The reason for this may be 
that the small number of blue stars detected toward the 
halo of NGC 253 are distributed over a large field. Wide 
field and deep imaging are needed to investigate in detail 
the star formation history in the halo of NGC 253, and its 
potential link with the observed galactic outflow. 

In the panel next to the color-magnitude diagram. Fig. 
3 shows the results of the edge detection flltering of the I- 
band luminosity function. The right panel of Fig. 3 shows 
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the logarithmic luminosity function. The power-law dis- 
tribution of the upper RGB stars is obvious, and has a 
clear break at Itrgb ~ 23.5 mag. Overplotted are the 
maximum-likelihood best fit model before and after con- 
volution with the photometric errors and application of 
incompleteness, shown as dashed and solid lines respec- 
tively. The bottom panel shows the posterior probabil- 
ity distribution of the TRGB magnitude derived from the 
Monte-Carlo Markov Chain analysis. Also shown are the 
68% and 95% confidence intervals. Both the Sobcl edge de- 
tector and the maximum likelihood model agree on the lo- 
cation of the I-band TRGB magnitude. For the maximum 
likelihood analysis, with a foreground reddening of E(B- 
V)=0.02 mag and fitting data in the range 22.8 < / < 24.2, 
we derive that the TRGB is located at /t_rgs — 23.5 mag. 
The Monte-Carlo Markov Chain analysis gives that the 
68% confidence interval spans the range 23.50 < I < 23.54 
and the 95% confidence interval is 23.48 < / < 23.79. The 
faint extention of the 95% confidence level is quite large, 
when the TRGB is well defined, as there are a variety of 
models with different amplitudes for the jump and slopes 
above and below the TRGB that are consistent with the 
data to within the 95% confidence interval. There is a 
family of solutions with a steep powerlaw at the bright 
end and a much smaller jump than the best-fit shown in 
the plot. The edge detection method locates the TRGB 
I-band magnitude at Itrgb = 23.51 mag with a boot- 
strap uncertainty of ±0.07mag. Following the procedure 
described in § 4, we find a Population II distance modu- 
lus (m - M)o = 27.59 ± 0.06(random)±0.16(systematic). 
This result is consistent with NGC 253 being a member 
of the nearby Sculptor group, located at the far side of 
the group; however, it is significantly larger than previ- 
ously reported values derived using the globular cluster 
luminosity function (Blecha 1986), and the brightest halo 
(Davidge & Pritchet 1990) and disk (Davidge, Le Fevre, 
& Clark 1991) stars. 

4.2. NGC 4244 

The edge-on galaxy NGC 4244 is a late-type (Scd) spi- 
ral; its total B-band magnitude is almost a magnitude be- 
low the peak of the field Sc luminosity function, making 
this galaxy less massive than a typical Sc. The radial light 
distribution is well represented by an exponential profile 
(Oiling 1996), with almost no dark lanes (i.e., internal ex- 
tinction is not important). 

The color-magnitude diagram for the observed halo 
field, shown in the left panel of Fig. 4, reveals a simple 
stellar population; the red giant branch is prominent and 
clearly dominates the halo light, and there seems to be 
no blue early-type stellar population in the halo of this 
galaxy. The red giant branch is relatively tight, indicating 
a tight metallicity distribution in the halo stars, with a 
clear tip at 7 « 23.9 mag. The stars above the tip are 
likely AGB stars. Also shown is the output of the edge 
detection filtering of the I-band luminosity function. 

The upper right panel of Fig. 4 shows the logarithmic 
luminosity function for all stars detected in the field. For 
NGC 4244, the TRGB appears as an inflection in slope, 
rather than a clear discontinuity. The power-law luminos- 
ity distribution of stars on the red giant branch is clear 
down to the magnitude where the effect of incompleteness 



starts to be important; the break of the luminosity func- 
tion due to the TRGB, at / « 23.9, is obvious. Again 
both the edge detection method and the maximum like- 
lihood model agree on the location of the TRGB I-band 
magnitude. Using the maximum likelihood analysis, and 
fitting data in the range 23.4 < I < 24.7, we find that 
the TRGB is located at Itrgb = 23.9 mag. The Monte- 
Carlo Markov Chain analysis gives a 68% confidence in- 
terval of 23.88 < / < 24.08 and a 95% confidence interval 
of 23.82 < / < 24.22. The edge detection method locates 
the TRGB at Itrgb = 23.88 mag. Resampling simula- 
tions give uncertainties of ±0.1 mag for the edge detection 
method. With a foreground extinction of E{B — V) = 
0.021 mag, we find the distance modulus to NGC 4244 to 
be (m - M)o = 27.88 ± 0.06(random)±0.16(systematic), 
corresponding to a heliocentric metric distance of D = 
3.82 Mpc. 

Our distance modulus is in agreement with the pre- 
viously published distance modulus for NGC 4244 of 
/z = 27.78 using the TuUy-Fisher relation (Aaronson et 
al. 1986), but significantly different from the distance de- 
rived by Karachentsev & Drozdovsky (1998), who used the 
mean B-band magnitude of blue supergiant candidates to 
derive a distance modulus of /U = 28.28. This method is af- 
fected by large uncertainties (Rozanski & Rowan- Robinson 
1994 and references therein). A preprint by Tikhonov & 
Galazutdinova (2005) analyzing a larger data set on this 
galaxy, including the field analysed in this paper, appeared 
after our paper was submitted. Their derived distance 
modulus using the tip of the red giant branch technique is 
28.16; which is consistent within our 2(t confidence interval 
but not our la confidence interval. The magnitude of the 
tip of the red giant branch derived by Tikhonov & Galazut- 
dinova (2005) is fainter than what is estimated here (see 
Fig. 3 of paper III of this series, and Fig. 4 of Tikhonov 
& Galazutdinova 2005). On the other hand, their esti- 
mated mean stellar metallicity is lower than our estimate 
(see paper II of this series). Tikhonov & Galazutdinova 
(2005) have not well documented their data reduction and 
the procedure they used to measure the distance modulus 
to NGC 4244 to trace back accurately the source of the 
observed discrepancy. 

4.3. NGC 4945 

NGC 4945 is an edge-on Sb galaxy, hosting a powerful 
nuclear starburst region, with a ring morphology (Moor- 
wood et al. 1996), powering a galactic superwind (Heck- 
man, Armus, & Miley 1990). Similar to NGC 5128, the 
prominent galaxy of the group to which NGC 4945 is be- 
lieved to belong, its optical image is marked by dust ex- 
tinction in the nuclear regions (Marconi et al. 2000; Lipari 
et al. 1997). 

The left panel of Fig. 5 shows the color-magnitude di- 
agram of the observed halo field stellar population. The 
red giant branch is wide and well populated, indicating a 
large spread of the halo stellar metallicities. Also shown 
is the output of the edge detection filtering of the I-band 
luminosity function. 

The right panel of Fig. 5 shows the logarithmic I-band 
luminosity function, compared to the maximum-likelihood 
best model derived by fitting the data in the range 
23 < / < 24.2. Again the upper red giant branch 
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power law distribution is evident with a break around 
I 23.5 — 23.6. The edge detection method and the 
maximum likeHhood model disagree slightly on the lo- 
cation of the I-band TRGB. Using the maximum likeli- 
hood technique to fit the stellar luminosity function in 
the range 23 < / < 24.2, we locate the TRGB mag- 
nitude at Itrgb = 23.63 mag, with a 68% confidence 
interval of 23.57 < I < 23.71 and a 95% confidence in- 
terval of 23.49 < / < 23.77. The edge detection method 
gives a brighter location of the TRGB at Itrgb = 23.48 
mag, with a bootstrap uncertainty of ±0.06mag. The 
average location of the TRGB is Itrgb = 23.55 mag. 
Using this value with a foreground reddening of E(B- 
V)=0.18mag, we find the distance modulus for this galaxy 
to be (m-M)o = 27.63±0.06(random)±0.16(systematic). 
This result is consistent with NGC 4945 being member of 
Centaurus A group, which has a centroid at a distance 
modulus of ^ = 27.81 ± 0.11 (3.66 ± O.lQMpc; Karachent- 
sev et al. 2002). 

4.4. NGC 4258 

NGC 4258 is a nearby barred spiral galaxy hosting an 
obscured active nucleus and a massive nuclear black hole 
inside a highly inclined thin gaseous disk (Miyoshi ct al. 
1995; Maoz 1995). Nuclear water maser sources were dis- 
covered moving in a Keplerian orbit around the central 
black hole of NGC 4258 (Watson & Walhn 1994). By re- 
solving these mascrs, measuring their acceleration (Green- 
hill et al. 1995), and monitoring their proper motions, 
Herrnstein et al. (1999) derived a distance modulus to the 
galaxy of (m - M)o = 29.28 ± 0.09. NGC 4258 is one of 
the only two galaxies with a geometrical distance estimate, 
the other galaxy being the Large Magellanic Cloud, with a 
geometrical distance measurement based on the light echo 
of SN 1987A (Panagia et al. 1991). Thus NGC 4258 is 
a unique object for a cross-check of extragalactic distance 
techniques. 

Recently, Newman et al. (2001) used a sample of 
Cepheids in NGC 4258 to report the distance modu- 
lus to this galaxy is (m - M)o = 29.47 ± 0.09 mag, 
with a distance modulus to the LMC of {m — M)o = 
18.50 mag. To do so, they used revised calibrations 
and methods for the HST Key Project on the Extra- 
galactic Distance Scale. This distance is 1.2 a larger 
the the maser determination. To reconcile the geomet- 
rical and the Cephcid distances to NGC 4258, they ar- 
gue that the true distance modulus to the LMC should 
be 18.31 ±0.11(random)±0.17(systematic), and should be 
used to calibrate the Cepheid brightness. Using the plan- 
etary nebulae luminosity function technique, Ciardullo et 
al. (2002) derived a distance modulus to NGC 4258 of 
(m - M)o = 29.42t°;^o. Similar to Newman et al (2001), 
they argue that the short distance scale to the LMC should 
be the correct one. On other hand, Caputo et al. (2002) 
argue that the theoretical metallicity correction, as sug- 
gested by pulsation models, is ~ 0.2 mag lower than the 
empirical correction adopted by Newman et al. (2001). 
Then if this theoretical correction is applied to account for 
the chemical abundance differences between the LMC and 
NGC 4258, it may reconcile the Cepheid and the maser 
distances to NGC 4258, with the long distance scale to 
the LMC, a scale supported by different empirical non- 



Cepheid methods (Panagia et al. 1991; Sakai et al. 2000; 
Cioni et al. 2000; Groenewegen & Salaris 2001). This de- 
bate raises again the well-known controversy between the 
short and long distance scales to the LMC. 

We take advantage of our halo field observations to esti- 
mate the distance modulus to NGC 4258 using the TRGB 
technique. Provided we restrict the ourselves toV — I < 2, 
this method is relatively insensitive to metallicity. The 
TRGB luminosity calibration does not depend on the as- 
sumed distance to the LMC, but rather is computed from 
theory and tested and calibrated via observations of globu- 
lar clusters in our galaxy. Because NGC 4258 is metal rich 
and the LMC is relatively metal poor, the TRGB distance 
to NGC 4258 is an important cross-check of the metallicity 
effect on the Cepheid-bascd extragalactic distance scale. 

The left panel of Fig. 6 shows the color-magnitude dia- 
gram of the observed halo field stellar population. Again 
the halo stellar population is dominated by red giant 
branch stars. The red giant branch is wide, indicating 
a large spread of the stellar metallicity in the halo. Also 
plotted is the output of the edge detection filtering of the I- 
band luminosity function. The right panel of Fig. 6 shows 
the logarithmic luminosity function for all stars observed 
in the field, compared to the maximum likelihood best fit 
model. The bottom panel shows the posterior probabil- 
ity distribution of the TRGB magnitude derived from the 
Monte-Carlo Markov Chain analysis. 

Again both techniques agree on the location of the 
I-band TRGB magnitude. For the maximum likeli- 
hood analysis, with a foreground foreground of E(B- 
V)=0.02mag and fitting data in the range 24.2 < I < 26, 
we locate the TRGB at Itrgb = 25.25 mag. The 
Monte-Carlo Markov Chain analysis gives that the 68% 
and 95% confidence intervals are 25.23 < I < 25.38 
and 25.19 < I < 25.46, respectively. The edge detec- 
tion technique locates the TRGB at Itrgb = 25.22 
mag with a bootstrap uncertainty of ±0.09 mag. Follow- 
ing the procedure described in § 4, we find a Population 
II distance modulus for this galaxy to be (m — M)o = 
29.32 ± 0.09 (random) ±0.1 5 (systematic), or a metric dis- 
tance of D=7.3 ± 0.3 Mpc. This is in superb agreement 
with both the maser and the metallicity-corrected Cepheid 
distances to the galaxy. Our distance modulus argue for 
the long distance scale to the LMC, where the true dis- 
tance modulus to the LMC is close to 18.50 mag. The 
results provide support for Hubble; constant Ho = 72 ± 8 
km s~^ Mpc~^, inferred by Freedman et al. (2001) and 
based on the long distance scale to the LMC, in contrast 
to other recent suggestions (Maoz et al. 1999; Ciardullo 
et al. 2002). 

5. CONCLUSION 

We report distance measurements using data obtained 
from a new program of deep HST imaging of the halos of 
nearby bright, edge-on spiral galaxies. 

We have used the TRGB as a distance indicator for four 
nearby galaxies, NGC 253, NGC 4244, NGC 4945, and 
NGC 4258, observed with HST/WFPC2. We employ two 
methods, edge detection and maximum likelihood fitting, 
to measure the TRGB I-band magnitude from the lumi- 
nosity function of RGB stars. Our results are in good 
agreement with prior distance measurements based on a 
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variety of methods, including Cepheid variables. However, 
we now have a common method of distance measurement 
for these galaxies. The uncertainty in the estimate of the 
luminosity function break is calculated using the bootstrap 
resampling technique, which makes no assumption about 
the underlying distributions. Using our TRGB distance 
to NGC 4258, in combination with the galaxy's geometric 
and new metallicity-corrected distances, we argue that the 
long distance scale to the Large Magellanic Cloud is cor- 
rect, and that the Frccdman et al. (2001) Hubble constant 
should not be revised. 

Finally, increasing the number of galaxies beyond the 
Local Group with accurate TRGB measurements will cer- 
tainly refine our understanding of the extragalactic dis- 



tance scale, and will help making the systematic compari- 
son between population I and population II distance indi- 
cators more meaningful. 

In future papers in this series, the same data set is used 
to derive abundance distributions in the halo stellar pop- 
ulations. 
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Table 1 

Distance modulus for the galaxy sample. Columns: (1) Galaxy name; (2) distance 

MODULUS estimates FROM THE LITERATURE; (3) THE METHOD USED TO ESTIMATE THE 

DIST_\.K('E.: (1) oris DISTANCE AIODFU'S ESTIMATE: (J)) AlISOLTTE I-DAM) AIACXITTDE . 



Galaxy 


(m - M)°" 


method 


(m - M)™'' 


Mi 


NGC 4258 


29.28 ± 0.06 


Maser ^ 


29.32 ±0.09 


-4.06 ± 0.04 




29.31 ±0.06 


Cepheids + Z-correction^ 








29.47 ±0.12 


Cepheids"^ 








29.421°:?^ 


PNLF* 






NGC 253 


27.30 


Globulars LF^'^ 


27.59 ±0.06 


-4.09 ±0.02 




< 26.8 


brightest halo and disk stars*'^ 






NGC 4244 


28.28 


brightest supergiants^° 


27.88 ±0.06 


-4.02 ±0.03 




27.78 


TuUy-Fisher relation^^ 






NGC 4945 


27.82 


Cen-A group centroid distance^^'^ 


27.63 ±0.06 


-4.06 ±0.02 



^Herrnstein et al. 1999 
^Caputo et al. 2002 
^Newman et al. 2001 
^Ciardullo et al 2002 
^Ferrarese et al. 2000 
•^Blecha 1986 
■^Puche & Carignan 1988 
^Davdige & Pritchet 1990 
^Davidge et al. 1991 
^^Karachentsev & Drozdovsky 1998 
^^Aaronson et al. 1986 
^^Karachentsev et al. 2002 
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Fig. 1. — Locations of the observed halo fields being studied in the paper. 




Fig. 2.— 



Evolution of photometric errors (in magnitude) as a function of the apparent I-band magnitude. 
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Fig. 3. — Left: color-magnitude diagram for stars observed in the NGC 253 halo field. The photometry has been corrected for foreground 
reddening. Center: the output of the weighted logarithmic edge detection filtering function. Right: Luminosity function on a logajrithmic 
scale. The data have been modeled with two power laws using a maximum-likelihood fit. Overplotted are the input best fit model (dashed 
line) and the model after convolution with the magnitude errors and application of incompleteness (solid line). The bottom panel shows the 
posterior probability distribution of the I-band TRGB magnitude derived from the Monte-Carlo Markov Chain analysis including the full 
range of the luminosity function free parameters. The extension of the 95% confidence toward fainter TRGB magnitudes interval represents 
a family of solutions that have steeper bright-end power-law slopes and smaller break-amplitudes than the best-fit model shown here. 
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Fig. 4.— Similar to Fig. 3 but for NGC 4244 
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V-I E(I) 23 23,5 24 24.5 

I(TRGB) 



Fig. 5.— Similar to Fig. 3 but for NGC 4945 
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Fig. 6.— Similar to Fig. 3 but for NGC 4258 



